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We have studied the v + V — > e + + e~ energy deposition rate in a rotating compact star. This 
reaction is important for the study of gamma ray bursts. The General Relativistic (GR) effects on 
the energy deposition rate have been incorporated. We find that the efficiency of the process is larger 
for a rotating star. The total energy deposition rate increases by more than an order of magnitude 
due to rotation. The dependence of this energy deposition rate on the deformation parameter of 
' the star has also been discussed. 
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I. INTRODUCTION 

Gamma Ray Burst (GRB) and its possible connection with neutrino production in compact stars is a field of high 
current interest. GRBs were first discovered in the late 1960s by U.S. military satellites Until recently, GRBs 
^*T* ' were perhaps the biggest mystery in high energy astronomy. To unravel this mystery, several satellite based detectors 
, have been employed for the observation of GRB. GRBs are separated in two classes; long duration bursts (long GRB) 
which last from 2 sec. to several minutes, with average duration of 30 seconds and short duration bursts (short GRB) 
with burst duration from few milliseconds to 2 seconds, average being 0.3 seconds @. 

Initial evidence that the long GRBs are associated with supernovae came from the study of GRB 980425 in 1998 
Q. This burst was tentatively linked to supernova SN 1998bw. The definitive proof came on March 29, 2003, when 
a relatively nearby burst, GRB 030329, produced an afterglow whose optical spectrum was nearly identical to a 
supernova Q. In contrast, there is a scarcity of information on short GRBs. Until recently, information was only 
available about the burst; the post-burst picture was not clear. Of late, the detection of afterglows in short GRBs and 
also precise localizations of different short GRBs have provided some more inputs for the study of these phenomena. 

On the otherhand, compact stars (neutron or quark stars) are objects formed in the aftermath of supernova. The 
central density of these stars can be as high as 10 times that of normal nuclear matter. Due to beta equilibration, 
a large number of neutrinos and antineutrinos may be produced inside the compact stars. These neutrinos and 
antincutrinos could annihilate and give rise to electron-positron pairs through the reaction vv — > e + e~. These e + e~ 
pairs may further give rise to gamma rays which could be a possible explanation of the observed GRBs. Hence it is 
very important to study the energy deposition in the vv annihilation process. 

Previous calculation of this reaction in the vicinity of a neutron star has been based on newtonian gravity 0, @] , 
i.e. (2GM/c 2 R) << 1, where M is the gravitational mass of the neutron star and R is the distance scale. The effect 
of gravity was incorporated in refs. 0, S|, but only for a static star. In our present calculation, we extend the basic 
premise to rotating stars. 

In this work first we will discuss the equation of state (EOS), the metric and the structure of the star. The effect 
of GR and also that of rotation will be discussed next and finally we will present our results. 
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II. STAR STRUCTURE 

The structure of the star is described by the metric given by Q 

ds 2 = -e 1+p di 2 + e 2a {dr 2 + r 2 d6 2 ) + e^r 2 sin 2 9{d<j, - udt) 2 (1) 

The four gravitational potentials, namely a, 7, p and cj are functions of and r only. Once these potentials are 
known, the observed properties of the star can be evaluted. The Einstein's equations for the three potentials 7, p 
and cj have been solved using Green's function technique [l(J EH • The fourth potential a has been determined from 
other potentials. All the potentials have been solved for both static as well as rotating stars using the 'rns' code; the 
details of the code may be obtained in ref. 

The solution of the Einstein's equations needs an EOS as an input. In the present work we have used a non-linear 
version of the Walecka model with TM1 [l3| parameter set. For comparison, we have also used a quark matter EOS 
obtained from standard Bag model with B 1 / 4 — 160MeV. These EOSs have been plotted in figure 1. The quark 
matter EOS is more stable at higher energy densities. Using these EOSs, we can get the structure of the star by 
solving the Einstein's equations. The solution of Einstein's equations also provide the density profiles of the star. 
These profiles, for a star rotating in the mass shedding limit i.e. with the Keplerian velocity, are plotted in fig. 2 and 
fig. 3. The neutrino annihilation process has been studied for both neutron and quark stars with the central energy 
density 1.2 x 10 15 gm / cm 3 . Fig. 2 shows the variation of energy density with the radius of the star for \ = cosO = 0, 
i.e. along the equator. Here 9 is the angle that radius vector makes with the polar axis. From fig 2. we can see that 
the energy density is high at the centre of the star, and as we move radially outwards the energy density decreases to 
a lower value. Fig. 3 shows the variation of the energy density with \- This is plotted for radius r = Z.bKm from the 
centre of the star. The curve shows that for X = 0, i.e. at the equator, the energy density is maximum and it falls 
off gradually as we approach the pole x = 1- We would like to mention at this stage that for the hadronic matter we 
have considered a thin crust. On the other hand, quark matter being self bound, no such crust has been used as for 
quark star. Hence the energy density at the surface of a quark star is much higher compared to that of a hadronic 
star (fig. 2). 
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FIG. 2: Variation of energy density along the equator of the star which is rotating with the Keplerian velocity. 
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FIG. 3: Variation of energy density with \ a t distance of 3.5Km from the centre of the star. 



III. GR EFFECTS 



As already mentioned above the energy deposition rate due to the process vv — > e + e in Newtonian gravity has 
been studied earlier 0, Q . The energy deposited per unit volume per unit time in Newtonian gravity can be given as 
[fill 
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(K r ) = / / L(Pv,r)fv(pi7, r)[a\v v - i^e^e-] ^ + £ " d 3 p u d 3 pi7, (2) 

where /„ and /u are the number densities of neutrinos and antincutrinos respectively in the phase space, v v is the 
neutrino velocity, and a is the rest frame cross section for the process vv — > e + e~. In eqn. [21 we have, 

DGj 

[a\v v -vjjle^} = — — (ei/e F -p„.pF) • (3) 

The above expression is Lorentz invariant. G 2 = 136 x lO~ 24 MeV~ 4 , and I? = 1 ± 4sm 2 6*vK + &sin 4 8w- Here 
sin 2 9w = 0.23. We will consider the plus sign only in the expression of D since this corresponds to the v e v e pairs. It 
has been assumed that the mass of the electron is negligible. Now we replace p v — e^fl^ and d 3 p L/ — t v d(. v d£l v , where 
f2„ is the unit direction vector and d^l v is the solid angle. So the rate of energy deposition is 

DG 2 r r 

= ^r e M / / + ep)44^ u dep (4) 

where the angular integration is represented by 

6(r) = / J (1 - Sl u Sy) 2 dn v dSlx?- (5) 

So the energy and angular dependences can be decoupled to make the evaluation simpler. 

Let us now consider the effects of GR on this process. The effect of GR on static star was studied by Salmonson 
and Wilson Q. In the present paper, we consider the static as well as the rotating star. The GR effect will modify 
both the energy and angular integrals and will provide us with a new lower limit of integration. 

Let us first consider the path of a zero mass particle, i.e. a null geodesic. For the metric considered here, one gets, 

, 1 dr~ 1 1 1 2u 11 

( ) 2 = 1 (Q) 

y r 2 e"t-Pd(t) b 2 e 2a {e"f+P - e"t~Pr 2 u> 2 ) b (e^+P - e^~Pr 2 uj 2 ) e 2a r 2 e a- P > w 

where r is the distance from the origin, <j> is the latitude and b is the impact parameter. From the above equation, 
one can immediately see that the geodesic explicitly depends on the gravitational potentials i.e. on the EOS and the 
frequency of rotation. If we now follow the technique used by p} and express this equation in terms of the angle 9 
between the particle trajectory and the tangent vector to the circular orbit, we get 

If the above expression is substituted in eqn. (6), we get a quadratic expression for b. 

b 2 - J^l 2 b~ 7n T \ =0. (8) 
e — v uj e p — t oj 

This equation can be solved to get 

ujr 2 ± re p 

b =—o 9 

e p — r oj 

b will be same for all the points lying on the same orbit. The above equation also implies a minimum photosphere 
radius, which we would denote by R, below which a massless particle (neutrino) emitted tangentially to the stellar 
surface (9r — 0) would be gravitationally bound. 

To carry out the angular integral in eqn. (5), we define A = sin9, and express the solid angle Q as 

ft = [A, (1 - A 2 )5 C os0, (1 - \ 2 )isin<j>] (10) 
dn = cos8d9dcj). (11) 



Using the same notation, we further obtain 
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where A„ and A p are for neutrinos and antineutrinos respectively, and 



* (r) = V e^-R^ R (13) 

Pr being the potential p at the photosphere. Carrying out the integration over X u and A-, we get 

2tt 2 

0(r) = —(1 - x(r)) 4 {x(r) 2 + Ax{r) + 5). (14) 

The gravitational effect modifies the temperature in the energy integral. The temperature T(r) in the energy integral 
eqn. (4) is the neutrino temperature at radius r. The temperature of the free streaming neutrinos at radius r in terms 
of their temperature at the minimum photosphere radius R can be deduced if we assume that the temperature varies 
linearly with redshift. Therefore 



In order to quantify the total e + e~ pair energy deposition, we define Q as the integral of q over proper volume. Q is 
a measure of the total amount of energy converted from the neutrinos to e + e~ pairs over the whole volume per unit 
time. So the total amount of local energy deposited via v + v — > e + + e~ reaction is given by 

r surface f 2a + (~j-p)/2 

Q = 2ir dr q{r)r 2 sinO d6. (16) 

Jr J VI -v 2 

with v = (fl — uj)rsin6e~ p , ft being the rotational velocity of the star and is the lorentz factor. 

Integration over the radial and angular variable gives us the total energy deposited per unit time by the vv pairs. 
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FIG. 4: Variation of the ratio of rate of energy deposition for general relativistic (GR) over Newtonian is plotted along the 
radial direction, for the hadronic star. 
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FIG. 5: Variation of the ratio of rate of energy deposition for general relativistic (GR) over Newtonian is plotted along the 
radial direction, for the quark star. 



trapping may increase the temperature considerably inside the star. The eqn. (|15p shows that the temperature of 
the star decreases as we move radially outwards because temperature, like energy, varies linearly with redshift. For 
the static star, we find that b is minimum for r — R = IKm for the hadronic star and r = R = \.2Km for quark 
star. For a rotating star due to high rotational velocity of the star, the shape of the star becomes oblate spheroid. 
The star is therefore no longer spherically symmetrical and it has an extra deformation parameter to describe its 
shape. We define this deformation parameter \ — cos8, along the vertical axis of the star. For the EOS, which we 
have taken the 'rns' code gives the rotational velocity of the star to be O = 0.62 x 10 4 s _1 for the hadronic star and 
Q = 1.4 x 10 4 s _1 for the quark star, the central energy density of both the stars being 1.2 x 10 15 gm/cm 3 . For such 
a configuration the star deforms to a oblate spheroid but the photosphere is found to be spherical. Photosphere is 
evaluated by minimising b for different x- The photosphere radius comes out to be r = R = 0.78 Km for the hadronic 
star and r = R = 0.6Km for the quark star. 

Figures 4 and 5 depicts the importance of GR on the energy deposition rate both for the static and the rotating 
compact stars. In both these figures we have divided the radius of the star into small bins of length 100m and 
integrated Q over those bins. The curves in fig. 4 show that for a static star, the energy deposition increases by a 
factor of 2 at the photosphere. As we go towards the surface the deposition increases and at the surface of the star 
this ratio becomes little over 3. For the rotating star, at the surface, this ratio is a little less than 8. For the hadronic 
star the two curves intersect. The ratio for the quark star is about 4 for the static case and reaches a much higher 
value of 20 at the surface when the star is rotating. From these two figures we conclude that the effect of rotation is 
to enhance the energy deposition by an order of magnitude both for the quark star and the hadronic star. Moreover, 
a comparison of the figures 4 and 5 also gives the effect EOS on the energy deposition rate. Since most of the other 
results are qualitatively similar, in the rest of the paper we will use hadronic star for further discussions. 

The figures 6 and 7 show the effect of GR and rotation on the energy deposition rate. In fig. 6 the variation of 
has been plotted as a function of \, f° r different r. Here, one can see that the ratio changes with the polar 
angle due to rotation. The deposition rate is somewhat pronounced at the pole compared to the equator, mainly near 
the surface of the star. In fig. 7 we have plotted the total energy deposition rate along the radial direction. The 
rate of total energy deposition is maximum at the photosphere and it decreases radially outwards, which is implied 
by the fact that the curves saturate to a maximum limit at higher radius. From this figure one can see that due to 
rotation the total energy deposition increases by an order of magnitude to 10 52 ergs/ s. It should be mentioned at 
this stage that the stars (both static and rotating) are constructed with equal central density (1.2 x 10 15 gm/cm 3 ). 
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FIG. 6: Variation of the ratio of rate of energy deposition with \i f° r three different radial points inside the neutron star. 
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FIG. 7: Variation of rate of energy deposition of general relativistic (GR) and Newtonian along the radius of the static and 
rotating neutron stars. 



mass sequences. The qualitative features remain the same and the quantitative results change only by a few percent. 
For comparison we should mention that Salmonson and Wilson 0, Q did the GR calculation for such a reaction on 
a static star using the Schwarzschild metric. They found that for type II supernova models the enhancement is by a 
factor of 2 at a rather larger radius R — 10M, increasing to a factor of 4 at R = 5M, where M being the total mass 
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argued that 0.5 x 10 53 ergs may be liberated from a star in thermal neutrinos within a few seconds |17j . 

To conclude, we have studied the effect of GR and rotation on the energy deposition rate of the reaction vv — ► e + e~. 
The effect of rotation on such a reaction has been considered for the first time. Our findings show that these effects 
enhance the energy deposition rate by more than an order of magnitude. This finding is very important in the 
context of describing the GRB associated with a compact, massive, rotating object. In our present calculation we 
have not considered any neutrino trapping. However the neutrinos may be trapped and raise the temperature of the 
star [3 EH) thereby further changing the energy deposition rate. Furthermore the consideration of non-equilabrated 
quark matter rather than equilibrated matter may enhance the energy deposition further. Such studies are in progress. 
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